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We discuss the role of the axial U(1)a symmetry in the chiral phase transition using the U {N/)n x 
U{Nj)l linear sigma model with two massless quark flavors. We expect that above a certain 
temperature the axial E/(1)a symmetry will be effectively restored as well as SU(Nf)n x SU(Nf)i,. 
Then we can construct a string-like static solution of the r\ string and a kink-like classical solution 
of the domain wall during the chiral phase transition. We give out the possible signals for detecting 
the domain wall in ultrarelativistic heavy-ion collisions. 
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Exploring the phase structure of quantum chromody- 
namics(QCD) is one of the primary goals of ultrarela- 
tivistic heavy- ion collisions. It is generally believed that 
at sufficiently high temperature there should be a transi- 
tion from ordinary hadronic matter to a chirally symmet- 
ric plasma of quark and gluons [l| . The order parameter 
for this phase transition is the quark-antiquark conden- 
sate. At temperature of about 150 MeV, Lattice QCD 
calculations indicate that this symmetry is restored 0- 
The order of the phase transition seems to depend on the 
mass of the non-strange u and d quarks, m u rs mj, and 
the mass of the strange quark m s , and at the tempera- 
ture on the order of 150 MeV, heavier quark flavors do 
not play an essential role. 

For Nf massless quark flavors, the QCD Lagrangian 
possesses a chiral U(Nj)r x U(Nj)l = SU(Nf)n x 
SU(N f ) L x U(l)v x U(1) A symmetry, here V = R + L, 
while A = R — L. However this symmetry does not ap- 
pear in the low energy particle spectrum, it is sponta- 
neously broken to the diagonal SU(Nf)v group of vec- 
tor transformation by a non-vanishing expectation value 
of the quark-antiquark condensate, (q R qL) ^ 0. This 
process involves Nj Goldstone bosons which dominate 
the low-energy dynamics of the theory. The U(l)v sym- 
metry is always respected and thus plays no role in 
the symmetry breaking pattern considered in the fol- 
lowing discussion. The axial U(1)a symmetry is broken 
to Z(Nf)A by a non- vanishing topological susceptibility 
[3 . Consequently, one of the Nf Goldstone bosons be- 
comes massive, leaving Nf — 1 Goldstone bosons. The 
SU(N f ) R x SU(N f ) L x U(1) A group is also explicitly 
broken by the effects of nonzero quark masses. 

As the temperature or the density of matter increases, 
it is expected that the instanton effects will rapidly dis- 
appear, the U{\)a symmetry may also be effectively re- 
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stored in addition to SU(Nf)n x SU(Nf)L- Since the 
chiral condensate (q R qL) ^ also breaks the U(1)a 
axial symmetry, there are only two possibilities: either 
the U(1)a symmetry is restored at a temperature much 
greater than the SU(Nf) R x SU(Nf)L symmetry or the 
two symmetries are restored at (approximately) the same 
temperature. Recent lattice gauge theory computations 
have demonstrated a rapid dropping of the topological 
susceptibility around the chiral phase transition, seem- 
ingly suggesting that the simultaneous restoration ex- 
ists 0, this is also supported by the random matrix 
models 5j. On the other hand, the fate of the U{1)a 
anomaly in nature is not completely clear since instan- 
ton liquid model calculations indicate that the topologi- 
cal susceptibility is essentially unchanged at T c jfj, also 
Lattice results obtained from the St/ (3) pure gauge the- 
ory show that the topological susceptibility is approxi- 
mately constant up to the critical temperature T c , it has 
a sharp decrease above the transition, but it remains to 
be different from zero up to ~ 1.2T C |7J. Additionally, 
other lattice computations which measure the chiral sus- 
ceptibility find that the U(1)a symmetry restoration is 
at or below the 15% level [| 0). 

Recently, the issue of finding signals for the restora- 
tion of chiral symmetry in ultrarelativistic heavy-ion col- 
lisions has received considerable attention. For exam- 
ple, the signals for the restoration of the SU (2) chiral 
symmetry associated with the a meson have been pro- 
posed in Refs.pd|p~l|]. In particular, signals for detecting 
the effective restoration of the U(1)a chiral symmetry in 
ultrarelativistic heavy-ion collisions have been invoked 
in[I3[I3El. 

On the other hand, in QCD, spontaneous symmetry 
breaking U(N f ) R xU(N f ) L -> U(N f ) v in the chiral limit 
allows for existence of topological string defects, the for- 
mation and evolution of such defects and their possible 
observable effects in ultrarelativistic heavy-ion collisions 
as well as in the early universe transition have been in- 
voked in Rcfs. 151 161 171. In this letter, we study the 
effects from effective restoration of the U(1)a symmetry 
by using the U(Nf) R x U(Nj)l linear sigma model with 
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two massless flavors. 

The Lagrangian of the U(Nf)n x [/(Nf)^ linear sigma 
model is given by |18| 



£($) = Tr(a A1 $ + 9 /i $-TO 2 $ + $)-Ai[Tr($+$) 



-A 2 Tr($+$) 2 +c 
+Tr[iJ($ + $+)]. 



:[det($) + det($ + )] 



(1) 



Where $ is a complex Nf x Nf matrix parametrizing the 
scalar and pseudoscalar mesons, 



$ = T a (t> a = T a (a a 



(2) 



with a a being the scalar (J p — + ) fields and 7r a being 
the pseudoscalar (J p = CP) fields. The Nf x Nf matrix 
H breaks the symmetry explicitly and is is chosen as 



H = T n h n 



(3) 



where h a are external fields, a = 0, 1, • • • , Nj — 1 and 
T a ,a ^ are a basis of generators for the SU(Nf) Lie 
algebra. T = 1 is the generator for the ^(1)^ Lie alge- 
bra. 

In the above model, the determinant terms correspond 
to the U(1)a anomaly, as shown by 't Hooft jjjj, they 
arise from instantons. These terms are invariant under 
SU(N f ) R x SU(Nf) L = SU(N f )v x SU(N f ) A , but break 
the U{1)a symmetry of the Lagrangian explicitly. The 
last term in Eq. Q which is due to nonzero quark masses 
breaks the axial and possibly the SU(Nf)y vector sym- 
metry explicitly. 

When h a = 0, c = 0, for m 2 < the global SU(N f ) v x 
U(Nf)A symmetry is broken to SU(Nf)v, and ($) de- 
velops a non-vanishing vacuum expectation value, (<J>) = 
To^o- Spontaneously breaking U(Nf)A beads to Nj 
Goldstone bosons which form a pseudoscalar, N 2 dimen- 
sional multiplet. However when h a = 0, and c ^ 0, the 
U(1)a is further broken to Z(Nf) by the axial anomaly, 
and SU(Nf) v x SU (Nj)a is still the symmetry of the La- 
grangian. A nonvanishing ($) spontaneously breaks the 
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1 di- 
mensional multiplet. The N 2 th pseudoscalar meson is no 
longer massless, because the U(1)a symmetry is already 
explicitly broken, e.g for Nf = 2, the rj meson is mas- 
sive compared to other pseudoscalar mesons. All these 
symmetry are in addition explicitly broken by non-zero 
quark masses making all the Goldstone bosons massive. 

In the present study, since we only concentrate on the 
effects of the effective restoration of the U{1)a symme- 
try, we can ignore the possible effects of the restoration of 
SU(2)ji x SU(2)l, this implies that we can forget it and 
a fields, keeping only the a and rj mesons which are usu- 
ally specified by the U(1)a phase. With this restriction 
on $, the effective Lagrangian we adopt here is 

£($) = Tr(9 A1 $+9 Al $-m 2 <I> + $)-Ai[Tr($ + $)] 2 

-A 2 Tr($+$) 2 + c[det($) + dct ($+)], (4) 



symmetry to SU{Nf)v, with the appearance of Nj 
Goldstone bosons which form a pseudoscalar, N 2 



where $ = ^(a + irj)\. In the following, we demonstrate 
that both a static string-like solution of the rj string and a 
static kink- like solution of the domain wall are expected 
to be produced during the chiral phase transition 1 . 

The rj string is a static configuration of the Lagrangian 
of Eq.(QJ with c = 0. In this case, during chiral symme- 
try breaking, the field a takes on a nonvanishing expecta- 
tion value, which breaks U{2)r x U(2)l down to f/(2)y. 
This results in a massive a and four massless Goldstone 
bosons. 

In our discussion of the rj string and domain walls it is 
convenient to define the new fields 



V2 



(5) 



The linear sigma model in Eq. Q with c = now can be 
rewritten as 



£ = (a^)*(a^)-A(0>--) 2 



(G) 



where v 2 



and A = Ai 



For static configura- 



tions, the energy functional corresponding to the above 
Lagrangian is 

E = J d 3 x[V0* + \(d>*4> ~ y ) 2 ], (7) 
and the time independent equation of motion is 

V 2 = 2A(0>- (8) 

The rj string solution extremising the energy functional 
of Eq.Q is given in Refs. [ ^ [ll f . 



(9) 



where p(r) = [1 — cxp(— fir)], the coordinates rand are 
polar coordinates in the x — y plane, the rj string is as- 
sumed to lie along the z axis and fi 2 = A^-tj 2 . The 
energy per unit length of the string is 



E = [0.7b + \og{fiR)]irv 2 



(10) 



For global symmetry in general the energy density of the 
string solution is logarithmically divergent, R is intro- 
duced as a cutoff which is taken to be O(fm) in the fol- 
lowing numerical calculation. 

In the case of c ^ 0, during chiral symmetry break- 
ing, the field a takes on a nonvanishing expectation 



1 For simplicity we consider here the configurations which are spec- 
ified by the U(1)a phase only. In considering non-abelian phases, 
there is another class of topological defects known as non-abelian 
strings llTl . the pion strings, which can also exist during the chiral 
phase transition ll5l . 
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value, which breaks SU(2) R x SU(2) L down to SU(2) V . 
This results in a massive a and three massless Goldstone 
bosons, in the same time the n meson is massive com- 
pared to other pseudoscalar mesons. Then the determi- 
nant term in Eq.Q can not be simplistically neglected 
during the chiral phase transition in nature, so that one 
of the appropriate description is no longer one of the n 
strings, but one of domain walls. Then in the following 
discussion we only consider the possible effects of domain 
walls and ignore the possible effects of the n string in the 
ultrarelativistic heavy-ion collisions. With the defination 
of new fields in Eq.lJSJ), the Lagrangian of Eq.QJ can be 
simplistically expressed as 

C = {d „(/))* (d"(f>) - m 2 4>*(t> + cRe(0 2 ) - A (</>*</>) 2 . (11) 

After defining c = am 2 , the potential takes the form 



V{<t>) = \((/>*<t)) 2 - m 2 (aRe((f> 2 ) - 0* 



(12) 



The limit a — > 00 corresponds to the maximum ex- 
plicit ^(l)^ symmetry breaking. In this limit, for re- 
alistic values of the a meson and the 7r meson mass(i.e., 
m 2 — c = constant), the 77 and ao mesons become in- 
finitely heavy and are thus removed from the spectrum of 
physics excitations, and U(2)r x U(2)l is identical to the 
0(4) model, there has no 77 strings and domain walls. For 
the chiral symmetry spontaneously breaking to occur, we 
always require a > 1 . In the following numerical calcula- 
tion, we take c = (386.79MeF) 2 , for other parameters we 
have Ai = -31.51, A 2 = 82.77 and to 2 = (263.83Me]£) 2 
corresponding to m a = AOQMeV and m„ = 547Mey|18j. 

For static configuration in Eq. Hllf) . the energy func- 
tional is given by 

E= I d 3 z[V0*V(/> + \{4>*4>) + m 2 (0*0) - am 2 Re(<j> 2 )}. 



(13) 

The corresponding equation of motion for the field <j) is 
V 2 + m 2 (a0*-^)-2A|0 2 |0 = O, (14) 
which accepts the static symmetric kink solution [l9|[2"oT| 



'(a-1) , r (a- 1) n 
<t m\ I — ^ — tanh[y mx\, (15) 

n = 0. (16) 
The thickness of this wall is approximately 

5 ~ (mVcT- "Tr 1 ~ 0.7 fm, (17) 
and the mass per unit area of the walls is 
2V2m 3 
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-(a - 1)5 ~ (129.273MeI/) d . (18) 



law is derived = e^ v d v (j), the associated charge of a 
configuration is N = J dxj° = <j>\ x =+oo — <f>\x=-oo, the 
presence of a kink with <f> in different vacuum at x — ±00, 
gives rise to a non-zero charge N and consequently indi- 
cates the stability of the configuration. Moreover, the 
form of the potential implies that the symmetric wall so- 
lution (within the domain wall the n — 0) is dynamically 
stable. We consider infinitesimal perturbations of the 
field n and check if the variation in the energy is positive 
or negative. Discarding terms of cubic and higher orders 
in ?7, we find 



E E(domainivall) ~t~ 



(19) 



where 



The stability becomes a consequence of a topological con- 
servation law. The topological current from which this 



SE = J d 3 x[^VnVn+ + l)m 2 n 2 + ^cr 2 n 2 }. (20) 

From the above equation, the term 8E in Eq. i|19|) is always 
positive, therefore, the domain walls of the Lagrangian 
Qis topologically stable and dynamically stable. 

In the Sine-Gorden model, the kink solutions are abso- 
lutely stable and such a stable domain wall will immedi- 
ately rule out by the cosmological constraint in general. 
In our case, the domain wall is only metastable in full 
theory since there are other dynamical fields correspond- 
ing to the remaining SU(2) generators (such as 7r and 
ao fields). However, one can show that these dynamical 
fields do not contribute to the domain wall background 
but simply remain in their vacuum states. Their fluctu- 
ations affect the overall energy density, but do not affect 
the properties of the domain wall such as the surface 
tension and so we can neglect their effects 21] . Then the 
domain wall can still be taken as classically stable object, 
and therefore, it decays through the quantum tunnelling 
process with exponentially large lifetime which is longer 
than any other time scales existing in the ultrarelativis- 
tic heavy ion collisions [22^. Then all the pions which are 
eventually emitted from such an object will be completely 
incoherent with the rest of pions. 

In the ultrarelativistic heavy-ion collisions, domain 
walls are expected to be produced durin g th e chiral phase 
transition. If a bubble wall is produced [22j, it exists for 
some lifetime and then decays into its underling fields, 
the a fields. We make the assumption that the size of 
the bubble wall should be around the size of the QGP 
formed at ultrarelativistic heavy-ion collisions. The ex- 
perimental observation of the domain wall bubbles can be 
carried out by using the Hanbury-Brown-Twiss (HBT) 
intensity interferometry of pionsE^] H3 • As pointed by 
Shuryak and Zhitnitsky in Ref.|22| if a bubble exists for 
enough long time(~5 fm) and then decays the bubble 
can be taken as an long-lived object. Therefore the pions 
from the bubble lead to the same effect of not produc- 
ing an HBT peak in two-pion spectra which is just as 
that of the long-lived hadronic resonances. To see this, 
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an effective intercept parameter, A, 
Bose-Einstein correlation function 2 



is introduced in 



C2(k ' K) = NMNM = * + V/ftW*,*), (21) 

where the effective intercept parameter and the correlator 
are given by 



<W/(p) = 



N c (p) 



1 2 



and 



R c {k,K) 



N c (p)+N h (p) 



S c (k,K) 



(22) 



S c (k = 0,K = p) 



(23) 



where k = (pi - pa) , K = (p a + pa)/2 , iV c (JV h ) is 
the one-particle invariant momentum distribution of the 
"core" ( and "halo" ) decayed pions respectively. S c is 
the Fourier transform of the one-boson emission function. 
The produced bubbles would given an additional factor 
j3 to the effective intercept. 



A' 



eff 
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PXeff, 



(24) 



where Ndomain-waii is the number of pions from the de- 
cay of domain wall bubbles. In RHIC Pb-Pb collisions if 
we take the radius of QGP phase as the domain wall bub- 
ble radius R ~ 6/jto 26J , then the domain wall bubble 
energy is about E domainwaU ~ 4irR 2 cu w 25GeV, . If all 
the energy accumulated in the wall will lead to the pro- 
duction of the a mesons (which will result in additional 
~ 60 mesons per event) one should expect a 40 7r + (or 
7r~) to be produced from the bubble wall in the central 
rapidity region. At RHIC energy the total number of pi- 
ons is about 1500, so the factor is about (3 ~ 0.85 . In the 
case of LHC Pb-Pb collisions the QGP radius is about 10 
fm|2^, this gives out j3 ~ 0.7—0.8. Thus we can use pion 
interferometry as a sensitive tool to detect this possible 
increase of the a production in ultrarelativistic heavy-ion 
collisions. 

In summary, we have discussed the possible effects of 
the restoration of the axial U(1)a symmetry during the 
chiral phase transition by using the U(Nf)n x U(Nf)r, 



linear sigma model with two massless quark flavors. It 
is emphasized that if the axial U{\)a symmetry is to be 
restored above the certain temperature, it is the domain 
wall rather than the t] string that is expected to be pro- 
duced and has a long lifetime then the time scale exist- 
ing in the ultrarelativistic heavy-ion collisions. These do- 
main walls will decay into the a mesons, and the increase 
of the a mesons can be viewed as a signal of restoration of 
the axial U{1)a symmetry in ultrarelativistic heavy-ion 
collisions. 

The authors wish to thank Michiyasu Nagasawa, 
Nicholas Petropoulos and Xinmin Zhang for useful dis- 
cussions and correspondence. 
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